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We report on a study of the decay B° —^ D'"^ujtv with the BABAR detector at the PEP-II B- 
factory at the Stanford Linear Accelerator Center. Based on a sample of 232 million BB decays, 
we measure the branching fraction B{B° ^ D*+lott~) = (2.88 ± 0.21(stat.) ± 0.31(syst.)) x 10"^. 
We study the invariant mass spectrum of the u>tt~ system in this decay. This spectrum is in 
good agreement with expectations based on factorization and the measured spectrum in r^ -^ 
uj-R~Vr. We also measure the polarization of the D*+ as a function of the uj-r~ mass. In the 
mass region 1.1 to 1.9 GeV we measure the fraction of longitudinal polarization of the D*"*" to be 
Fl/F = 0.654 ± 0.042(stat.) ± 0.016(syst.). This is in agreement with the expectations from heavy- 
quark effective theory and factorization assuming that the decay proceeds as B" — > D*+/9(1450), 
p(1450) -^uj-K-. 

PACS numbers: 13.25.Hw, 12.39.St, 14.40.Nd 



I. INTRODUCTION 

Factorization is a powerful tool to describe hadronic 
decays of the B-meson. According to factorization, the 
matrix element of four-quark operators can be written as 
the product of matrix elements of two two-quark opera- 
tors Q- Thus, the process b — > cW* , W* — > qq' (where 
q = d or s, q' — u or c) can be "broken up" into two 
pieces, the b —* c transition and the hadronization from 
W* ^ qq' decay. 

Ligeti, Luke, and Wise have proposed an elegant test of 
factorization 2J. In this test, data from r —^ Xu, where 
AT is a hadronic system, is used to predict the properties 
of S — > D*X (see Fig. ^. If A" is composed of two or 
more particles not dominated by a single narrow reso- 
nance, factorization can be tested in different kinematic 
regions. 






U/ X = UJTT 



^r B'>rv 



is large (~ 0.3%), this decay provides a good laboratory 
for the study of factorization. 

The branching fraction for i3° -^ D*^ujTr^ has been 
measured by the CLEO collaboration, using a sample of 
9.7 million BB pairs collected at the T(4S') resonance, 
to be (2.9 ± 0.3(stat.) ± 0.4(syst.)) x lO^^ 4]. They also 
extracted the spectrum of rrij^, the square of the invari- 
ant mass of the ujtt system. This spectrum is found to 
be in agreement with theoretical expectations U ■ In ad- 
dition, the CLEO collaboration studied the related de- 
cay B -^ Dlutt and concluded that this decay is domi- 
nated by the broad p(1450) intermediate resonance; i.e., 
B -^ Dp{14:50), p(1450) -^ ujtt. Assuming that this inter- 



mediate state also dominates in 5° ~> D*'^u 



factor- 



ization can be used to predict the polarization of the D* 
with the aid of heavy-quark effective theory (HQET) and 
data from semileptonic i?-decays J3- These predictions 
are in agreement with the CLEO result for the longitu- 
dinal polarization fraction, Tl/T = (63 ± 9)% 4]. 

In this paper we study the decay i?° -^ D*'^ujtt~ with 
a larger sample of B decays than available in the original 
CLEO study. We present measurements of the branching 
fraction, the Dalitz plot distribution, the m^ spectrum, 
the mu'T: distribution, and the D* polarization as a func- 
tion of nix ■ 



FIG. 1: Feynman diagrams for r ^ lj-kv and B^ -^ D*^ui-k . 

In the event that AT is a multi-body system, it is pos- 
sible that some fraction of the hadronic system could be 
emitted in association with the B — > £)'*■* transition in- 
stead of the hadronization from W* — > qq' decay. In the 
case of A" = wtt", the pion must come from the W* to 
conserve charge. It is unlikely that the uj could be pro- 
duced from the lower vertex in Fig. ^0) 13- Furthermore 
the LOTi^ state is not associated with any narrow reso- 
nance, so that a wide range in wtt" invariant mass can 
be studied. As the branching fraction for B^ -^ D*^ujtt^ 
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II. THE BABAR DATASET AND DETECTOR 

The results presented in this paper are based on 232 x 
10^ T(4S') -^ BB decays, corresponding to an integrated 
luminosity of 211 fb~^. The data were collected between 
1999 and 2004 with the BABAR detector Q at the PEP-II 
B Factory at SLAG. In addition a 22 fb~^ off-resonance 
data sample, with center-of-mass energy 40 MeV below 
the T{AS) resonance, is used to study backgrounds from 
continuum events, e~^e~ -^ qq {q = u, d, s, or c). 

Charged-particle tracking is provided by a five-layer 
double-sided silicon vertex tracker (SVT) and a 40-layer 
drift chamber (DCH), operating within a 1.5-T magnetic 
field. Energy depositions are measured with a CsI(Tl) 
electromagnetic calorimeter (EMC). Charged particles 
are identified from ionization energy loss (dS/dx) mea- 
surements in the SVT and DCH, and from the observed 
pattern of Cherenkov light in an internally reflecting ring 



imaging detector. 



III. ANALYSIS STRATEGY 

Starting from the set of reconstructed charged tracks 
and energy deposits within the EMC, we select events 
that are kinematically consistent with i?° — > D*^ujTr~ in 
the foUowing decay modes: D*^ -^ I?°7r+, with D'^ -^ 
K~Tr^, A'^TT+TT+TT^, Or K^tt^tt'^, and uj —^ tt+tt^tt'^. 
Charge-conjugate modes are implied throughout this pa- 
per. 

In the reconstruction chain, the invariant mass require- 
ment on the 7r+7r~7r° system that forms the uj candidate 
is kept loose. We then select "signal" or "sideband" can- 
didates depending on whether the reconstructed tt^tt^tt'^ 
mass is consistent with the w hypothesis. Kinematic dis- 
tributions of interest, such as the mj^ spectrum, are ob- 
tained by subtracting, with appropriate weights, the dis- 
tributions for signal and sideband events. This subtrac- 
tion accounts for all sources of backgrounds, including 



D*+TT-TT+TT-TT°, 



on a statis- 



backgrounds from B 

tical basis. This is because, as we will demonstrate in 

Section IVII background sources with real uj decays are 

neghgible. 

The event reconstruction efficiency is determined from 
simulated Monte Carlo events, where the response of the 
BABAR detector is modeled using the GEANT4 p] pro- 
gram. Efficiency-corrected kinematic distributions are 
obtained by assigning a weight to each event. This weight 
is equal to the inverse of the efficiency to reconstruct 
that particular event given its kinematic properties. This 
procedure, which is independent of assumptions on the 



dynamics of the B^ 
Section rvTrl 



D 



*+, 



decay, is discussed in 



IV. EVENT SELECTION CRITERIA 

The event selection criteria are optimized based on 
studies of off-resonance data, and simulated BB and con- 
tinuum events. 

Photon candidates are constructed from calorimeter 
clusters with lateral profiles consistent with photon show- 
ers and with energies above 30 MeV. Neutral pion can- 
didates are formed from pairs of photon candidates with 
invariant mass between 115 and 150 MeV and energy 
above 200 MeV, where the 7r° mass resolution is 6.5 MeV. 
In order to improve resolution, tt'^ — > 77 candidates are 
constrained to the world average 7r° mass |3l- 

The kaon-candidate track used to reconstruct the Z?" 
meson must satisfy a set of kaon identification criteria. 
The kaon identification efficiency depends on momentum 
and polar angle, and is typically about 93%. These re- 
quirements provide a rejection factor of order 10 against 
pions. For each D" -^ iir~7r+7r° candidate, we calculate 
the square of the decay amplitude (|^P) based on the 



kinematics of the decay products and the known prop- 
erties of the Dalitz plot for this decay [9|. We retain 
candidates if |j4p is greater than 2% of its maximum 
possible value. The efficiency of this requirement is 91%. 
Finally, the measured invariant mass of D^ candidates 
must be within 15 MeV of the world average D° mass [a| 
for 1)0 -^ K-TT+ andD° ^ X-tt+tt-tt-, and 25 MeV for 
Z)" -^ K^tt^tt'^. The experimental resolution is about 6 
MeV for D" ^ K-tt+, K-n+n-n-, and 10 MeV for 

We select D*+ candidates by combining Z?° candidates 
with an additional track, assumed to correspond to a 
pion. We require the measured mass difference Am = 
m{D*+) - TOpO) to be between 143.4 and 147.4 MeV. 
The resolution on this quantity is 0.3 MeV with non- 
Gaussian behavior due to the reconstruction of the low 
momentum pion from D* decay. 

In the rest frame of the 5°, as m^ increases the D*"*" is 
produced with decreasing energy. At high mj^ , or equiva- 
lently low Z3*+ energy, the reconstruction efficiency drops 
as cos ^D ^ 1 , where 9d is the angle between the daugh- 
ter D° and the direction opposite the flight of the B^ in 
the D*~^ rest frame. We exclude the region of low accep- 
tance (cos6'd > 0.8 for 8 < m^ < 9 GeV^ cos 6Id > 0.6 
for 9 < m^ < 10 GeV^, and cos 61c > 0.4 for mj^ > 10 
GeV^) from our event selection. The effect on the final 
results is very small, as will be discussed in Section IVlIII 

We form lu candidates from a pair of oppositely- 
charged tracks, assumed to be a tt+tt^ pair, and a tt*^ 
candidate. In order to keep signal and sideband can- 
didates (see Section IIII|I we impose only the very loose 
requirement that the invariant mass of the to candidate 
be within 70 MeV of the world average oj mass 8] . (The 
natural width of the uj resonance is 8.5 MeV and the 
experimental resolution is 5.6 MeV.) 

In order to reduce combinatoric backgrounds, we im- 
pose a requirement on the kinematics of the uj decay [lOJ . 
This is done by first defining two Dalitz plot coordinates: 
X = 3To/Q - 1 and y = \/3(T+ - T-)/Q, where T±,o 
are the kinetic energies of the pions in the uj rest frame 
and Q = T+ + T- + Tq. Next, we define the normal- 
ized square of the distance from the center of the Dalitz 
plot, i?2 = (X2 + Y'^)/{X^ + Y^^), where Xb and Yj, are 
the coordinates of the intersection between the kinematic 
boundary of the Dalitz plot and a line passing through 
(0,0) and {X,Y). Since the Dalitz plot density for real 
UJ decays peaks at i? = 0, we impose the requirement 
R < 0.85. This requirement is 93% efficient for signal 
and rejects 25% of the combinatorial background. 

We reconstruct a _B-meson candidate by combining a 
D*~^ candidate, an uj candidate, and an additional neg- 
atively charged track. A B-candidate is characterized 
kinematically by the energy-substituted mass ttt-es = 

1 /(is + po ■ pbY / E"^ — p^g, where E and p denote energy 

and momentum measured in the lab frame, the subscripts 
and B refer to the initial T(4S') and B candidate, re- 
spectively, and s represents the square of the energy of 
the e+e" center of mass (CM) system. For signal events 



we expect rriES ~ ^b within the experimental resolu- 
tion of about 3 MeV, where Mb is the world average 
B mass |a|. In the same fashion, the energy difference 



/^E 



Eh 



^Y^, where the asterisk denotes the CM 
frame, is expected to be nearly zero for signal B decays. 

The A£' resolution is approximately 25 MeV in the 
K^TT^TT^ mode and 20 MeV in the other modes, with 
non-Gaussian tails towards negative values due to energy 
leakage in the EMC. We select B candidates with a D^ -^ 
K~TT^TT^ as long as — 70 < AE < 40 MeV, and we require 
-50 < A£' < 35 MeV for the other modes. 

In order to further reduce the number of events from 
continuum backgrounds we make two additional require- 
ments. First, we require |cos0s| < 0.9, where 9b is the 
decay angle of the B candidate with respect to the e~ 
beam direction in the CM frame. For real B candidates, 
cos 9b follows al — x^ distribution, while the distribution 
is essentially flat for B candidates formed from random 
combinations of tracks. Second, we impose a require- 
ment on a Fisher discriminant [ij designed to differenti- 
ate between spherical BB events and jet-like continuum 
events. This discriminant is constructed from the quanti- 
ties Lq = J2iPi ^-iid -^2 = J2iPi cos^ a*. Here, p^ is the 
magnitude of the momentum and a^ is the angle with 
respect to the thrust axis of the B candidate of tracks 
and clusters not used to reconstruct the B, all in the CM 
frame. The requirements on | cos 9b \ and the Fisher dis- 
criminant are 95% efficient for signal and reject nearly 
40% of the continuum background. 



The reconstruction of the B^ 



D 



*+, 



decay is im- 



proved by refitting the momenta of the decay products 
of the B^ , taking into account kinematic and geomet- 
ric constraints. The kinematic constraints are based on 
the fact that their decay products must originate from 
a common point in space. The entire decay chain is fit 
simultaneously in order to account for any correlations 
between intermediate particles. 

If more than one B candidate is found in a given event 
with Toes > 5.2 GeV, and passes selection requirements, 
we retain the best candidate based on a x^ algorithm 
that uses the measured values, world average values, and 
resolutions of the Z?° mass and the mass difference Am. 
We omit the u candidate mass information from arbitra- 
tion in order to avoid introducing a bias in the lo mass 
distribution, since this distribution is used extensively 
throughout the analysis. 



V. EVENT YIELD 

In Fig. 121 we show the toes distribution for candidates 
with reconstructed 7r^7r~7r'' mass (rriuj) in the signal and 



,PDGl 



<20 



sideband regions, which are defined as |to^ 

MeV and 35 < \m^ - m^^°\ < 70 MeV respectively, 

where m^^'^ is the world average m mass Q • 

The TTiES distribution for the rrii^ signal region has been 
fitted to the sum of a threshold background function [IJ 
and a Gaussian distribution centered at Mb- The dis- 
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FIG. 2: rriES distributions for candidates with reconstructed uj 
mass in the signal (points) and sideband (shaded histogram) 
regions. The distribution for events in the sideband region 
has been rescaled to match the expected background in the 
rriuj signal region. The fitted function is described in the text. 



tribution for the rriaj sideband region demonstrates the 
presence of a background component, which peaks in 
Toes but not in to^;, that is not well described by the 
threshold function. Monte Carlo studies indicate that 
approximately one-third of this component is due to sig- 
nal events where the uj is mis-reconstructed. These are, 
for example, events where one of the pion tracks in the 
UJ decay is lost and is replaced by a track from the de- 
cay of the other B in the event. The remaining two- 
thirds of the Toes peaking background component is due 
to i?° — > D*+7r+7r^7r°7r^ events. 

We extract the event yield from a binned x^ fit of the 
m^ distribution for events with toes > 5-27 GeV. The 
data distribution is modelled as the sum of a Voigtian 
function and a linear background function. (The Voigtian 
is the convolution of a Breit-Wigner with a Gaussian res- 
olution function.) The width of the Breit-Wigner is fixed 
at 8.5 MeV, the world average width of the ui. The mass 
of the LU, the Gaussian resolution term, and the parame- 
ters of the linear function are free in the fit. 

The rrii^ distribution and the associated fit is shown 
in Fig. 121 The yield, defined as the number of events in 
the Voigtian with |to^ - m^'^'^\ < 20 MeV, is 1799 ± 
87 events. The Gaussian resolution returned by the fit 
as well as the mean of the Breit-Wigner are consistent 
with the value we find in Monte Carlo simulations of 
B''* -^ D*^ojTT~ events. In Fig. E|we also include the m^ 
distribution for events with 5.20 < toes < 5.25 GeV 
(the Toes sideband). This background distribution has 
been scaled to the number of background events expected 
from a fit to the ttt-es distribution where we require \m^ — 
ml^^^l < 70 MeV. The difference between the number of 
observed events away from the ?n„ peak and the number 
of background events predicted from the toes sideband 
is due to the background component that peaks in toes- 

The validity of the yield extraction relies on the as- 
sumption that the background is linear in to^^, and, most 
importantly, that there are no sources of combinatoric 
backgrounds that include real cu decays. The results 
shown in Fig. |31 imply that there is no significant com- 



> 400 

<U 

S 350 

CM 

;~-^300 

iS250 

a|200 
UJ 



100 
50 



: 




A 




r 




/\ 


'kiliJll^miliit^A+i 


^ 


* J**^ 




^Mf^m^ 


^r^^' 




nn-"^ 











0.72 0.74 0.76 0.78 



0.8 



0.82 0.84 

m,(GeV) 



FIG. 3; Distribution of reconstructed rric^ for events with 
JTiES > 5.27 GeV (points) and events with 5.20 < rriES < 5.25 
GeV (shaded histogram). The superimposed fit is described 
in the text. The events from the rriES sideband have been 
scaled to the expected background from an rriES fit to events 



with \m^ 
figure) . 



< 70 MeV (i.e., the range shown in this 



ponent of real lo decays in the background. To verify 
this, we have examined and fit the rriuj distribution for 
data events in the tties sideband as well as the distribu- 
tion for Monte Carlo simulations of BB events, excluding 
B^ -^ D*^LUTT^. Wc find that the distributions are well 
modelled by linear functions. There is no evidence of a 
real u) component in the background. We estimate that 
this component can affect the yield extraction of Fig. |3| 
at most at the few percent level. 

We also divide our dataset into three independent sub- 
datasets, according to the three D° decay modes that we 
consider. The fits to these sub-datasets yield consistent 
results. 



VI. BACKGROUND SUBTRACTION 

In this work we are interested in studying a number 
of kinematic distributions for S" -^ D*'^lljtt~, such as 
the m^ distribution, where mx is the invariant mass of 
the iOTT system. The measurements of these distributions 
need to account for the presence of background in the 
sample and for the fact that the signal reconstruction 
efficiency is not constant over the Dalitz plot. 

We use distributions for w sideband events to remove 
the effects of the background in the cu signal region on 
a statistical basis, and we use Monte Carlo simulations 
to correct for efficiency effects. This is accomplished as 
follows: 



1. The simulation of B' 



■0 



D 



*+, 



events is used to 



calculate the signal reconstruction efficiency e{x), 
where x is the set of quantities that specify the 



kinematics of a given event. The procedure used to 
determine e(x) is discussed in Section lVfll 

2. In the absence of background, we would calculate 
the number of events corrected for efficiency in a 
given bin of mj^ as 



N{mj,) 



E 

signal 



e{xi) ' 



(1) 



where the sum is over signal events in a given rnj^ 
bin and Xi is the set of kinematic quantities for the 
i-th event in the sum. 

3. As mentioned above, the background subtraction is 
performed using the rn^ sideband. Thus, Eq. ^is 
modified to be 



Nimj.) 



1 



Signal 



E 



1 



7 ^ — ' e[x 

sideband 



(2) 



where the first sum is just as before, while the sec- 
ond sum is over w-mass sideband events in the given 
bin of rnj^ and Xj represents the set of kinematic 
quantities for the j-th event in the sideband event 
sample. The same efficiency is used for both the 
signal and sideband event samples. The factor of 
I is needed to adjust for the relative size of the lj 
signal and sideband regions. The additional factor 
of /3 is ideally equal to one, and it is introduced 
to correct for any possible bias in the background 
subtraction procedure, as will be discussed below. 

The allowed kinematic limits for some variables, such 
as rn^, are not the same for lu signal and sideband events. 
Therefore, the values of these variables for events in the 
UJ sideband region are linearly rescaled so that their kine- 
matic limits match the kinematic limits for events in the 
LU signal region. This procedure is necessary to avoid 
the introduction of artificial structures in background- 
subtracted distributions for these variables near the kine- 
matic limits. 

Wc test the sideband subtraction algorithm on a 
number of background samples such as Monte Carlo 
BB events and data events in sidebands of wes and 
AE. These tests are performed using the efficiency 
parametrization discussed in Section IVIII We find 
that background-subtracted kinematic distributions in 
the background samples show no significant structure. 
One sample distribution is shown in Fig. 0] We find a 
small bias in the extraction of the background-subtracted 
yields if the parameter (3 in Eq.|21is set to unity. To cor- 
rect for this bias we set /? = 0.975, with an estimated 
systematic uncertainty of ±0.0f 0. 
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FIG. 4: (a) Efficiency-corrected m^ distributions for events 
from the niES sideband with reconstructed rriui in the signal 
and sideband regions (arbitrary units). The distribution for 
events in the sideband region has been scaled by a factor of | . 
(b) Background subtracted rux distribution for events from 
the tues sideband (arbitrary units). This distribution has 
been obtained by subtracting the two distributions in (a). 



VII. EFFICIENCY PARAMETRIZATION 

The process of interest (S" -^ D*^ujtt^) is the three- 
body decay of a pseudoscalar particle into two vector 
particles and a pseudoscalar particle. We parametrize the 
reconstruction efficiency as a function of five variables: 

1. d, an index that labels the decay mode of the I?°; 
i.e., D^ — > K^TT^, A'^7r+7r+7r~, or K^tt^tt'-'; 

2. Ei^, the energy of the u! in the 5° rest frame; 

3. Ed', the energy of the D* in the 5° rest frame; 

4. cos 9d, the cosine of the decay angle of the D*; 
i.e., the angle between the D^ and the direction 
opposite the flight of the S" in the £**+ rest frame. 

5. cos a, the cosine of the angle between the vector 
normal to the uj decay plane and the direction op- 
posite the flight of the B^, measured in the u; rest 
frame. 

Note that two other variables would be needed to fully 
describe the kinematics of the decay chain. These are 
the angles that define, in addition to cos6'_d and cos a, 
the orientation between the decay planes of the D* and 
the u! and the decay plane of the i?°. Monte Carlo studies 
show that the reconstruction efficiency is independent of 
these two variables. The E^ and Ed-^ variables are the 
usual Dalitz variables used to describe three-body decays. 
Because of energy-momentum conservation the E^^ and 
Eu» variables are equivalent in information content to 
the squared invariant masses of the D*tt (rnj^..,^) and ojtt 
(rnj^) systems respectively. 



The efficiency is then parametrized as 

e{xi) = e{E^, Ed', cos 9 D,\ cos a\;d) (3) 

= e'{E^,ED';d) •ci(£'^,|cosa|) ■ C2{Ed' , cos 6 D;d). 

The functions e', ci, and C2 are extracted from Monte 
Carlo simulations and tabulated as a set of two dimen- 
sional histograms. As an example, the e' distribution for 
events with D'^ -^ K^i:^ is given in Fig.|Sl 
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FIG. 5: The e'(_D" ^ K-ti+, E^, Ed-) distribution for B" 
D*^u}-K~ Monte Carlo events. 



The efficiency parametrization is validated using sam- 
ples of Monte Carlo signal events. These samples are 
generated with a variety of ad-hoc kinematic properties; 
e.g., different polarizations for the D* and the uj, different 
shapes of the rn?^ distribution. In all cases we find that 
the shapes of kinematic distributions are well reproduced 
after the efficiency correction. 

We use the following method to estimate the effect of 
the finite statistics of the Monte Carlo sample. We gen- 
erate a set of 400 new e', ci, and C2 templates based on 
the nominal templates obtained from Monte Carlo sig- 
nal events. If the measured efficiency in a given bin of 
the nominal template is /i ± cr, the corresponding effi- 
ciencies in the new templates are drawn from a Gaus- 
sian distribution of mean /i and standard deviation cr. 
Then, the measurement of any quantity of interest {e.g.. 
m 



X 



) is repeated 400 times, according to Eq. El using 
the new templates. The spread in the results obtained 
from events reconstructed in data is a measure of the 
systematic uncertainty due to the finite number of avail- 
able Monte Carlo events. This spread is then added in 
quadrature to the statistical uncertainty of our results. 



We observe a small bias in the total number of recon- 
structed signal events obtained from the efficiency cor- 
rection. This is due to the fact that, although the un- 
certainty on e{xi) is Gaussian, the factor l/e{xi) used in 
the efficiency correction procedure (Eqs. ^ ^ does not 
obey Gaussian statistics. As a result, after applying the 
efficiency correction, the total number of reconstructed 
events tends to slightly overestimate the true value. 

In order to quantify this bias on the nominal result due 
to the finite number of Monte Carlo signal events, we first 
determine the mean of the total number of reconstructed 
signal events in data for the 400 new efficiency templates. 
This mean differs from the nominal result by a few per- 
cent (S). We then repeat the procedure described above 
using events reconstructed from signal Monte Carlo. We 
use the results of these Monte Carlo studies to describe 
the bias as a function of S. We find that after applying 
the efficiency correction and subtracting the m^ side- 
band, the total number of events reconstructed using sig- 
nal Monte Carlo exceeds the true value by (0.6 ± 0.4) • 5. 
We correct our final results by this amount. 



VIII. RESULTS 

We use the procedure outlined above, with one addi- 
tional correction, to extract the branching fraction, the 
m^ distribution, the Dalitz plot distribution, the mw-n 
distribution, and the polarization of the D* as a func- 
tion of mx- The one additional correction accounts for 
the region of phase space with low acceptance that was 
excluded from the analysis. This region corresponds to 
values of cos 9d near 1 for low E'^ , or equivalently high 
m\ . This correction factor varies between approximately 
1.2 at TO^ = 8 GeV^ and 1.6 at m\ = 11 GeV^. Since 
most of the data is at rrix < 4 GeV^, the combined effect 
of this correction is quite small; it amounts to an increase 
of less than 1% relative to the measured branching frac- 
tion. 

For the branching fraction, we find B{B^ -^ 



D*+LU7r-) = (2.88±0.21(stat.)±0.31(syst.))xl0"^. The 
total systematic uncertainty of 10.8% arises from the fol- 
lowing sources: 

• The uncertainties in the branching fractions of the 
D*, D, and ui: 5%. 

• The uncertainty in the reconstruction efficiency of 
neutral pions at BABAR, which is estimated to be 
3% per tt". This amounts to a 6% uncertainty for 
events reconstructed with D^ -^ K~Tr~^Tr^, and 3% 
for the other modes. Combining these modes, the 
systematic uncertainty from this source is 4.3%. 

• The uncertainty in the reconstruction efficiency for 
charged tracks. From a variety of control samples, 
this is estimated to be 0.6% (0.8%) for each track of 
transverse momentum above (below) 200 MeV. In- 
cluding the uncertainty for the low momentum pion 



produced in the D*~^ decay, we obtain a systematic 
uncertainty of 5.3%. 

• The uncertainty in the efficiency of the kaon parti- 
cle identification requirements: 2%. 

• The uncertainty due to the limited Monte Carlo 
sample size in the efficiency calculation: 3.8%. 

• The uncertainty on the quantity (3 in Eq. 12 2.6%. 

• The uncertainty in the efficiency of the various 
event selection criteria, estimated to be 4.3%. 

• The uncertainty in the number of BB events in the 
BABAR event sample: 1.1%. 

• The uncertainty in the correction due to the re- 
moval of events at high cos 9d and small Eu* : 
0.3%. 

Some of these systematic uncertainties vary as a func- 



tion of 



For example, the uncertainty on the cor- 



rection due to removing a region of {Ed*,cos9d) phase 
space is only relevant to events with rnj^ above 8 GeV^. 
A portion of the systematic uncertainty due to limited 
Monte Carlo sample size also varies as a function of m|- . 
Therefore, quantities measured as a function of m^ in- 
clude a common scale uncertainty of 10.5% and a sys- 
tematic uncertainty that varies with m^ and is typically 
below a few percent. 

The nij^ distribution, normalized to the semileptonic 
width r{B ~> D*£i^) |i|, is shown in Fig. A scale un- 
certainty on our result of 11.3% is not shown. This uncer- 
tainty combines a 4.2% uncertainty in r{B -^ D*tv) with 
the 10.5% uncertainty from the sources listed above. The 
bulk of the data is concentrated in a broad peak around 
m\ « 2 GeV^, in the region of /9(1450) -^ uir. Our distri- 
bution agrees well in both shape and normalization with 
predictions based on factorization and r decay data J13l | 
in the region m\ < 2.8 GeV^ covered by the tau data. 

The background-subtracted and efficiency-corrected 
Dalitz plot is shown in Fig. [7| One notable feature of 
the decay distribution is an enhancement for D*tt masses 
near 2.5 GeV [mjj,^ ^ 6.3 GeV^). The enhancement 
occurs in the region where one expects to find a broad 
J = 1 D** resonance (D'l) that decays via S-wave to 
£'*7r. Thus, this enhancement could be due to the color- 
suppressed decay B^ — > D'iLj, followed by D'^ — > £>*+7r~. 

In Fig. |S1 we show the background-subtracted and 
efficiency-corrected D*n mass distribution for events 
away from the p(1450) peak, fitted to the sum of a 
fourth order polynomial and a relativistic Breit-Wigner. 
In this figure, in order to remove the contribution from 
the /9(1450), we have required cos 9d* < 0.5, where 0£). is 
the angle between the momentum of the D* in the D*tt 
rest frame, and the fiight direction of D*tt system. We 
use the cos6'£i* variable rather than m^ to remove the 
p(1450) contribution because the distribution in cos^d* 
is uniform for S-wave D[ -^ D*tt decay. The yield of 
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FIG. 6: (a) Data rrix distribution normalized to the semilep- 
tonic width r{B -^ D*lv). The inner error bars reflect the 
statistical uncertainties on the data. The total error bars in- 
clude the m^-dependent systematic uncertainties. A common 
11.3% scale systematic uncertainty is not shown, (b) Same as 
(a) but zoomed-in on the low nXx region, where comparisons 
based on factorization and r data can be made. Also shown 
here are the results from the CLEO analysis [j. 



possible B^ -^ D[uj events in Fig. |H| can then be easily 
extrapolated to the full kinematic range. Furthermore, 
by subdividing the dataset in bins of cos Od* we can test 
the S-wave decay hypothesis. 

The fitted mass and width of the Breit-Wigner in Fig.|Sl 
are 2477±28 MeV and 266±97 MeV, respectively These 
values are consistent with the parameters of the broad 
D'l -^ D*7r resonance observed by the Belle collaboration 
in B ^ D'^TT decaySj^ m = 2427 ± 36 MeV and T = 
3841^^^ ± 74 MeV fi3|. We have also spht the data set 
of Fig. |S1 into three equal-sized bins of cos^d*. We find 
that the fitted amplitude of the Breit-Wigner component 
is the same, within statistical uncertainties, in the three 
data sets. This is consistent with expectations for an 
S-wave D[ -^ D*it decay. 

If we assume that the enhancement for D*7r masses 
near 2.5 GeV is actually due to i3° -^ D'^ui, D[ -^ 
D*~^'!T~ , we extract the branching fraction 



B{B' 



dO 



-^ D[lj) X B{D[ -^ D*+TT 
(4.1 ±1.2 ±0.4 ±1.0) 



X 10" 



(4) 



In this measurement, the first uncertainty is statisti- 
cal, the second uncertainty is from the uncertainties in 
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FIG. 7; Background-subtracted and efficiency-corrected 
Dalitz plot for B — > D*uj-k. The relative box sizes indicate the 
population of the bins. Black boxes indicate positive values, 
white boxes indicate negative values, which can occur because 
of statistical fluctuations in the subtraction procedure. 
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FIG. 8: Background-subtracted and efficiency-corrected D*-k 
mass distribution with cos 9d' < 0.5. The superimposed fit 
is described in the text. 



the last uncertainty arises from the uncertainties on the 
choice of the nonresonant shape in Fig. [S] (10%) and 
the uncertainties in the parameters of the D[ resonance 
(22%). This branching fraction has been obtained from 
fitting the sample of events with cos6'£i' < 0.5, and scal- 
ing up the result by a factor of ^ . This procedure neglects 
interference effects between i3° -^ D[uj and B'^ -^ D*ujtt. 
The branching fraction in Eq. 0] is comparable to the 
branching fractions for B^ — *■ D^*'>^uj Q. Also, we see 
no evidence for decays into the two narrow D** reso- 
nances at 2420 and 2460 MeV. This is in contrast to the 



measurement, and color-favored i? — > D**"7r decays, where the three _D 
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modes contribute with comparable strengths, and where 
the B^ —f D[t:^ branching fraction is one order of mag- 
nitude smaUer than that of B^ -^ D^*'^^tt^ . 

The presence of B^ ^ D[uj would affect the compari- 
son of the data with the theoretical predictions of Fig.El 
As can be seen in Fig. [T] S" -^ D[llj would mostly con- 
tribute at high value of rrv^, while the factorization test 
can be carried out only where the t data is available; 
i.e., for rnj^ < 3 GeV'^. Based on the estimated branch- 
ing fraction of B^ -^ D'lUj, and neglecting interference 
effects, the contribution of B^ -^ D[uj to the m^ < 3 
GeV^ distribution would be less than 5%. 



If the decay B° 



D 



*+, 



proceeds dominantly 



through B° -^ D*+p{U50), p(1450) -^ ujtt- , a measure- 
ment of the polarization of the D* can provide a further 
test of factorization and HQET 15] . The angular distri- 
bution in the D*+ -^ Z?°7r+ decay can be written as a 
function of three complex amplitudes Hq (longitudinal), 
and i?-|_ and -ff_ (transverse), as 



dV 



d cos 9d 



(xA\Ha\^cos^eD- 



|i/+|2 + |i/_|2)sin2 0z,, (5) 



where 9d is the decay angle of the D* defined above. 
The longitudinal polarization fraction r^/r, given by 



Tl 

r 



^4 



H. 



(6) 



can then be extracted using Eq.Elfrom a fit to the angular 
distribution in the decay of the D* . 

We divide our dataset in ranges of m^, and per- 
form binned chi-squared fits to the efficiency-corrected, 
background-subtracted, I?*-decay angular distributions. 
In these measurements, nearly all of the systematic un- 
certainties discussed above cancel. As a result, the m?x- 
dependent uncertainty due to the finite Monte Carlo sam- 
ple is the dominant systematic uncertainty, and typically 
results in an uncertainty on Tl/T at the few percent level. 
We also include a systematic uncertainty due to the pa- 
rameter /3 in Eq. |21 This uncertainty is about one order 
of magnitude smaller. 

The measured longitudinal polarization fractions as a 
function of mx are presented in Table I. Near the mean 
of the /9(1450) resonance (1.1 < mx < 1.9 GeV), we find 
Tl/T = 0.654±0.042(stat.)±0.016(syst.). This result is 
in agreement with the previous result in the same mass 
range from the CLEO collaboration, Tl/F = 0.63 ±0.09. 
It is also in agreement with predictions based on HQET, 
factorization, and the measurement of semileptonic B- 
decay form factors, Fl/F = 0.684 ± 0.009 Hg, assum- 
ing that the decay proceeds via B^ -^ _D*+p(1450), 
/9(1450) -^ LUTT~ . These results are shown in Fig.^ 



IX. CONCLUSIONS 



We have studied the decay i3° -^ D*'^u)tt with a 
larger data sample than previously available. We mea- 



TABLE I: Results of the D* polarization measurement in bins 
of mx- The first uncertainty is statistical and the second is 
systematic. 



mx range 


(GeV) 


Fl/F 


below 1.1 




0.458 ± 0.189 ± 0.059 


1.1 - 1.35 




0.779 ± 0.062 ± 0.020 


1.35 - 1.55 




0.733 ± 0.071 ± 0.024 


1.55 - 1.9 




0.435 ± 0.102 ± 0.040 


1.9 - 2.83 




0.656 ± 0.182 ± 0.077 
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FIG. 9: The fraction of longitudinal polarization as a func- 
tion of m?x, where X is a vector meson. Shown (as a tri- 
angle) is the B'^ -^ D'^uj-K^ polarization measurement for 
events with 1.1 < mx < 1.9 GeV {mx = rrip,, where 
p' = p(1450)), as well as earlier measurements (indicated by 
open circles) of B° -^ D*+ p- ^, B° ^ D*+ D*' 0, and 
B — > D*'^ D*~ [l3 ). The shaded region represents the predic- 
tion (± one standard deviation) based on factorization and 
HQET, extrapolated from the semileptonic S" — > D*'^l~v 
form factor results [l9(|. 
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sure the branching fraction to be B{B^ -^ Z3*+ 
(2.88±0.21(stat.) ±0.32(syst.)) x lO'^. 

The invariant mass spectrum of the lott system is found 
to be in agreement with theoretical expectations based on 
factorization and t decay data. The Dalitz plot for this 
mode is very non-uniform, with most of the rate at low 
LOTT mass. We also find an enhancement for D*7r masses 
broadly distributed around 2.5 GeV. This enhancement 
could be due to color-suppressed decays into the broad 
D[ resonance, B^ -^ D[uj, followed by D[ -^ D*+n-, 
with a branching fraction comparable to B^ -^ D^*^^uj. 

We also measure the fraction of D* longitudinal po- 
larization in this decay. In the region of lott mass be- 
tween 1.1 and 1.9 GeV, where one expects contribu- 
tions from S" -^ D*+p{U50), p(1450) -> ujtt~ , we find 
Fl/F = 0.654 ± 0.042(stat.) ± 0.016(syst.), in agreement 
with predictions based on HQET, factorization, and the 
measurement of semileptonic -B-decay form factors. 
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